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Lattice dynamics and electron-phonon coupling in transition metal diborides 
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The phonon density-of-states of transition metal diborides TMB2 with TM = Ti, V, Ta, Nb and 
Y has been measured using the technique of inelastic neutron scattering. The experimental data are 
compared with ab initio density functional calculations whereby an excellent agreement is registered. 
The calculations thus can be used to obtain electron-phonon spectral functions within the isotropic 
limit. A comparison to similar data for MgB2 and AIB2 which were subject of prior publications as 
well as parameters important for the superconducting properties are part of the discussion. 
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The discovery of superconductivity in MgB2 at the 
unexpected high temperature of 39 K has renewed the 
interest in diboride compounds. The well known AIB2 
structure (P6/mmm) is formed with quite a number of 
transition metals (TM) 0. However, a clear proof of su- 
perconductivity could not be established for any of these 
candidates. Obviously the influence of lattice faults on 
the transport properties is important. A unique feature 
for MgB2 is the depopulation of B 2p — a bands due to 
a particular Mg-B interaction Q]. It has been shown by 
band structure calculations for the TM-diborides that Ep 
moves up into the TM 4d-states with the consequence 
that these compounds behave more like normal metals 
0- Most important for the transport properties is the 
electron-phonon interaction in these systems. Very re- 
cently results from point contact spectroscopy have been 
reported for TM = Zr, Nb, and Ta. The authors conclude 
on an only modest electron-phonon interaction parame- 
ter for these samples A general investigation of the 
phonon spectra of TM-diborides has so far not be per- 
formed. In the present paper we confront measurements 
of the generalized phonon density of states for selected 
compounds with TM = Ta, V, Nb, Ti, and Y with ab 
initio density functional theory (DFT) calculations. Re- 
sults for electron-phonon spectral functions are derived 
and placed in context to previous investigations of super- 
conducting MgB 2 . 

All diborides were prepared from stoichiometric mix- 
tures of TM elements and amorphous 11 B. We have cho- 
sen the less absorbing n B in view of the intended in- 
elastic neutron scattering (INS) experiments. The mixed 
powders were compacted to pellets and arc melted. The 
final polycrystalline samples showed metallic brightness 
and were found to be single phase by x-ray diffraction. 
All of our samples failed to show superconductivity above 
4K although transition temperatures for the pure ele- 
ments of Nb and Ta with 9.2 K and 4.39 K are known. 
Our INS experiments were performed on the IN6 time- 
of-flight spectrometer at the HFR in Grenoble, France, 
with an incident neutron energy of 4.75 meV at 300 K in 
the upscattering mode. A high chopper speed of 201 rms 
and focusing in the inelastic region was used to improve 
the resolution. The generalized phonon density of states 



TABLE I: Structural and scattering parameters for selected 
diborides. Distances are given in A. The values in brackets 
apply to optimized geometries found in the DFT calculations 
after energy optimization. M-B and B-B denote the shortest 
metal-boron and boron-boron distances, respectively, and m 
is the TM mass number, a /m is the neutron scattering cross- 
section over the mass in barn/amu. 

MgB 2 A1B 2 TaB 2 VB 2 NbB 2 TiB 2 YB 2 

— a 1084 31)09 3l)8 2^998 3^09 31)38 3^290" 

(3.056) (2.965) (3.08) (2.979) (3.093) (2.998) (3.254) 

c 3.522 3.262 3.265 3.056 3.3 3.23 3.835 

(3.622) (3.232) (3.272) (2.995) (3.337) (3.188) (3.830) 

M-B 2.504 2.383 2.421 2.297 2.477 2.409 2.716 

B-B 1.781 1.737 1.799 1.772 1.794 1.773 1.914 

m 24.312 26.982 180.948 59.42 92.906 47.90 88.91 

a/m 0.151 0.056 0.033 0.098 0.067 0.085 0.087 



(GDOS) has been calculated from the recorded intensi- 
ties integrated over a scattering region from 14° to 114°. 
For the data evaluation we have applied multi-phonon 
corrections in a self-consistent procedure. The GDOS 
implies a weighting of vibrational modes of by a/m (scat- 
tering cross-section over the mass, see Tab.Pl. For com- 
parison the scattering power of 11 B is 0.525 barn/amu. 

Fig. 1 shows a comparison of GDOS spectra for all 
investigated TM-diborides. Each spectrum is normal- 
ized to 1. The most significant common feature is a gap 
around 40 to 50meV. When analyzing the respective in- 
tensities we find that the high-frequency band is essen- 
tially composed of boron vibrations. The gap is thus a 
consequence of the large mass difference, which leads to 
a decoupling of transition metal and boron vibrations. A 
clearly different behavior is observed for MgB 2 and A1B 2 
where due to the smaller mass of the metal atoms a hy- 
bridization of metal and boron modes occurs. In the low- 
frequency band of the transition metal diboride spectra 
we observe distinct shifts which however do not simply 
follow a y/m relation. This indicates that there are appre- 
ciable differences in the B-TM coupling. This conclusion 
is confirmed by the concomitant variations of interatomic 
distances (Tab.QJ. Somewhat more surprising is the fact 
that significant frequency shifts are equally registered for 
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FIG. 1: The experimental generalized phonon density-of- 
states for various transition metal diborides. AIB2 and su- 
perconducting MgB2 are shown for comparison. Very signif- 
icant changes are related to structural differences. Signature 
of strong electron-phonon coupling is the dramatic broaden- 
ing of the highest energy B peak visible for MgB2 but not for 
the other candidates. 



the peak at the high-frequency end of the spectra where 
the lowest frequency is found in YB2. The F,2 g and 
Bi 5 modes, which are important for superconductivity 
in MgB2, lie within this region where neighboring boron 
atoms move against each other. Strong electron-phonon 
coupling is expected to cause a significant broadening. 
This broadening is visible in the spectrum of MgB2 but 
not in those of the other compounds. The low frequency 
of the peak in YB2 alone can, therefore, not be taken as 
a signature of strong electron-phonon coupling. A sim- 
pler explanation is offered by comparatively weak B-B- 
interactions. As can be seen from Tab.|]]the B-B distance 
is exceptionally large for this compound. If any, then the 
G(uj) for NbB 2 shows signatures of a stronger electron- 
phonon coupling, a compound for which some authors 
report superconductivity. 

First principals density-functional calculations give pa- 
rameter free insight into the electronic structure of all of 
these compounds. They in addition yield the lattice dy- 
namical properties. Corresponding results are shown in 
Fig. 2. For the calculations we used the mixed basis pseu- 
dopotential method, which is described in some detail in 
Ref. 0, 0. Norm-conserving pseudopotentials for V, 
Ta, Nb, Y and Ti were constructed according to the de- 
scription of Bachelet-Hamann-Schluter jjjla], whereas for 
boron a Vanderbilt-type potential was created @- The 
mixed basis scheme uses a combination of local functions 
together with plane waves for a representation of the va- 
lence states. The good agreement between calculated and 
experimental structural parameters (see Tab. P) as well 
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FIG. 2: Calculated GDOS spectra for all diborides shown in 
Fig. 1. The proper a/m values have been applied for a direct 
comparison to the experimental spectra. 



TABLE II: Calculated V point frequencies in meV. 





MgB 2 


AIB2 


TaB 2 


VB 2 


NbB 2 


T1B2 


YB 2 




40.5 


36.6 


52.6 


60.6 


52.0 


65.5 


44.1 


A 2 „ 


50.2 


52.1 


61.9 


62.1 


60.5 


66.4 


45.3 


E 2s 


70.8 


125.0 


100.6 


114.9 


98.4 


112.8 


75.5 


Big 


87.0 


61.3 


68.8 


69.6 


69.8 


70.0 


75.1 



as between calculated and experimental phonon GDOS 
(Figs. 1 and 2) proves the validity of our description. 

A very unique feature of electron-phonon coupling in 
MgB2 is the downshift of the in-plane E2 9 mode well be- 
low the out-of-plane Bi g mode. Such an inversion of the 
usual sequence of mode frequencies is not encountered in 
any of the other diborides as shown in Tab-ITTl where we 
compile the calculated T-point frequencies of the boron 
vibrations. 

The present calculations prove capable of describing 
the lattice dynamics of various diborides with quite dif- 
ferent electronic transport properties. For the TM di- 
borides, the isotropic character of the d-bands has the 
consequence that these compounds behave electronically 
very much like normal 3D metals Q ■ Although the band 
structures show strong similarities, the Fermi level Ep 
moves up and down within the region of the transition 
metal d-bands. This variable band filling leads to ap- 
preciable differences in the various interatomic distances. 
The E2 9 frequencies exhibit an expected inverse correla- 
tion with the in-plane lattice constant, which indicates 
that the B-B covalent bond is electronically rather simi- 
lar for all these compounds. On the contrast, the out-of- 
plane boron mode Bi g is very insensitive to changes in 



3 




ENERGY [meV] 

FIG. 3: Isotropic Eliashberg functions for the investigated 
diborides as obtained in the present calculation. 



the lattice constants. 

Somewhat exceptional in our series of transition metal 
diborides is only YB2. Peculiarities in the electronic 
structure lead to a drastic increase of both the a and 
c lattice constants. As a consequence of the increased B- 
B bond length, the longitudinal force constant of the B-B 
bond is strongly reduced, and becomes even smaller than 
the transverse coupling. This is the reason for the ob- 
served softening of the high-frequency B-modes in YB2 . 

In order to get a more quantitative picture of the 
electron-phonon coupling within the series of diborides 
we have calculated the isotropic Eliashberg function 
a 2 ¥(uo) which describes the phonon mediated pairing in- 
teraction (Fig. 3) . Details of the computations are out- 
lined in Ref. pg. Again MgB2 proves to be singular. The 
extremely strong coupling at high-phonon frequencies is 
known to be connected to the E2 9 boron in-plane mode. 
It is found as a general trend that the region of strong 
coupling for TM diborides shifts down to lower energies 
as found for conventional superconductors thus leaving 
hope only for modest T c 's. The existence of 2D- and 3D- 
type Fermi sheets derived mainly from boron a and 7r 
bands, respectively, has been shown to be important for 
a quantitative understanding of the transport properties 
of MgB 2 [13 • The 2D Fermi surface is connected to the 
hole doping of boron 2p bands and does not seem to be 
important for transition metal diborides. 

We have also calculated the relevant parameters for 
the description of the superconducting properties within 
the Eliashberg formalism. The results are compiled in 
Tab. II I II A common feature unfavorable for supercon- 
ductivity found in all transition metal diborides is the 
comparatively low LOi og . This quantity represents the ef- 
fective average frequency of the coupling modes and sets 



TABLE III: Average coupling constant A, effective average 
phonon frequency uii og (in meV), and estimates of T c from the 
linearized gap equation. N(0) denotes the density of states at 
the Fermi energy (per unit cell and spin). 





N(0) 


A 


<+>log 


T c (/x*=0.13) 


MgB 2 


0.335 


0.73 


60.9 


21.7 


A1B 2 


0.184 


0.43 


49.9 


2.3 


TaB 2 


0.452 


0.79 


25.8 


10.6 


VB 2 


0.592 


0.28 


44.1 


< 1 


NbB 2 


0.520 


0.67 


30.5 


8.4 


TiB 2 


0.179 


0.10 


52.9 




YB 2 


0.560 


0.46 


37.4 


2.4 



the energy scale for the pairing interaction. Its small 
value as compared to MgB2 indicates that the pairing 
interaction is mainly mediated by the TM vibrations and 
not by the boron modes. The isotropic coupling constant 
A shows a variation from weak to intermediate coupling 
strength. In Tab. IIIII we have also included values for 
the superconducting transition temperature as obtained 
by solving the linearized isotropic gap equations. The 
T c values depend strongly on the screening properties 
expressed by n* and drop appreciably for \i* > 0. The 
values shown for a typical metallic screening of /1* =0.13 
thus should only be taken as indicators of the general 
trend. Furthermore, as mentioned before, for MgB2 it is 
necessary to go beyond the isotropic limit and take into 
account the multigap structure introduced by the par- 
ticular Fermi surface geometry for a proper quantitative 
description of the pairing state. Anisotropy is expected 
to play a minor role for the coupling in the low-frequency 
region of the transition metal vibrations that are domi- 
nated by the 3D-mctal bands. 

The values of A and 0Ji og are most favorable for MgB2- 
The next promising candidates according to the present 
calculations are TaB2 and NbB2, for which some authors 
find superconductivity at 9.5 K and 6.4 K, respectively 
Pj . Differences in stoichiometry might be able to explain 
different experimental results. As far as a comparison 
to recent point-contact measurements is possible we can 
state that the maxima in our G(u>) correspond well to the 
peaks found in the second derivative of the I — V char- 
acteristic. There is also agreement in a strongly reduced 
electron-phonon coupling strength for transition metal 
diborides with respect to MgB2. We, however, want to 
mention that the reported very small numerical values for 
Xpc should not be compared directly to our values since 
point contact experiments do not measure the McMillan 
A. 

To summarize, we have presented measurements of the 
generalized phonon density-of-states for various transi- 
tion metal diborides. The observed trends in the lattice 
dynamics can be well ascribed by first principles calcu- 
lations and do not exhibit indications of strong electron- 
phonon interaction. This view is supported by theoreti- 
cal calculations of the electron-phonon coupling. While 
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they do not exclude the possibility of superconductivity the outstanding properties of MgB2 among the class of 
with low T c mediated by TM vibrations, they underline diborides. 
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